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Abstract 

Some issues of inverting asteroseismic frequency data 
are discussed, including the use of model calibration and 
linearized inversion. An illustrative inversion of artificial 
data for solar-type stars, using least-squares fitting of a 
small set of basis functions, is presented. A few details of 
kernel construction are also given. 
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1. Introduction 

The accumulated experience in helioseismology of invert- 
ing mode frequency data provides a good starting point 
for asteroseismic inversion. In some rather superficial ways 
the circumstances of the two may seem very different: 
helioseismologists can use many more mode frequencies 
than will ever be possible for a more distant star, so that 
the resolution that can be achieved in helioseismic inver- 
sion is beyond the grasp of asteroseismology (an exception 
may be situations such as where modes are trapped in a 
narrow range of depth within the star) . Another difference 
is that global parameters of the Sun such as its mass, ra- 
dius and age are much better known than they are for 
other stars; hence the structure of the Sun is constrained 
a priori much more accurately than it is for a distant star, 
even if the input physics had been known precisely, which 
of course is not the case for the Sun or for other stars. 

In more fundamental ways, though, helioseismic and 
asteroseismic inversion are much more similar than they 
are different. Helioseismology has not always been blessed 
with such a wealth of data, and helioseismologists in the 
early days of their subject learned the value of model 
calibration and asymptotic description (in particular of 
low-degree modes) for making inferences about the Sun. 
They also learned some of the dangers and limitations 
of drawing inferences from real data. The modal proper- 
ties in many asteroseismic targets will be similar to the 
low-degree solar modes. Moreover, the principles of inver- 
sion are the same in both fields: obtaining localized infor- 
mation about the unseen stellar interior; assessing resolu- 
tion; taking account of the effects of data errors; assessing 
what information the data really contain about the ob- 



ject of study; judiciously adding additional constraints or 
assumptions in making inferences from the data. 

This paper touches on a few of these points with regard 
to inverting asteroseismic data, including the usefulness 
of optimally localized averaging (OLA) kernels and model 
calibration using large and small frequency separations for 
solar-type oscillation spectra. Much more detail on those 
two topics can be found in the papers by Basu et al. (2001) 
and Monteiro et al. (2001), both in these proceedings. 

The inversion problem predicates that one is able to 
perform first the forward problem. In the present con- 
text, that means the computation of the oscillation data 
(the mode frequencies) from the assumed structure of the 
star. Solving the forward problem always involves approx- 
imations or simplifying assumptions, because we cannot 
model the full complexities of a real star. The inverse pro- 
cedure, inferring the structure (or dynamics) of the star 
from the observables, is ill-posed because, given one solu- 
tion, there will almost invariably formally be an infinite 
number of solutions that fit the data equally well. The art 
of inverse theory is in no small part concerned with how to 
select from that infinity of possibilities. Actually, in vari- 
ous asteroseismic applications to date, e.g. to r] Boo and 
certain 6 Scuti stars, the problem is rather that one has 
so far been unable to find even a single solution that fits 
the data (see Christensen-Dalsgaard, Bedding & Kjeldsen 
1995, Pamyatnykh et al. 1998). This may indicate that the 
some of the approximations made in the forward problem 
are inappropriate; it could also indicate that the errors 
in the data have been assessed incorrectly. In the present 
work, unless otherwise stated, we assume both that the 
forward problem can be solved correctly and that the sta- 
tistical properties of the data errors are correctly known. 

2. Model calibration 

A relatively straightforward approach to the inverse prob- 
lem is model calibration. At its conceptually simplest, this 
entails computing the observables for a set of models, pos- 
sibly a sequence in which one or more physical parameters 
vary through a range of values, and choosing from among 
that set the one that best fit the data. 'Best' here is of- 
ten taken to mean that model which minimizes the chi- 
squared value of a least-squares fit to the data. Given to- 
day's fast computers, searches through large sets of mod- 
els are feasible, either blindly or by using some search 
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algorithm such as genetic algorithm or Monte Carlo. Ap- 
proaches similar to this have been undertaken for white 
dwarf pulsators (Metcalfe, Nather & Winget 2000; Met- 
calfe 2001) S Scuti stars (Pamyatnykh et al. 1998), and 
pulsating sdB stars (Charpinet 2001). Depending on the 
way in which the method is applied, the resulting model 
may be constrained to be a member of the discrete set of 
calibration models, or could lie "between" them if inter- 
polation in the set of models is permitted. 

Model calibration is powerful and dangerous. It is pow- 
erful because it allows one to incorporate prejudice into 
the search for a solution and, suitably formulated, it can 
always find a best fit. It is dangerous for the same reaons: 
perhaps unwittingly on the part of the practitioner, it 
builds prejudice into the space of solutions that is con- 
sidered. Also, even a satisfactory fit to the data does not 
mean that the solution model is necessarily like the real 
star: one can make a one-parameter model calibration to 
a single datum, but it is unclear what aspects of the re- 
sulting stellar model the datum is actually able to con- 
strain. Examples of prejudice that model calibration may 
incorporate, for good or ill, are the choice of physics used 
to construct the models, and possible assumptions about 
the smoothness of the structure of the star. An insidi- 
ous problem is that, if the approximations in the forward 
modelling introduce errors into the computed model ob- 
servables, this can introduce a systematic error into the 
result of the model calibration. The Sun and solar-type 
stars provide a good example. Here the near-surface struc- 
ture and the mode physics in that region are poorly mod- 
elled at present, and the simple approximations made in- 
troduce a systematic shift in the computed frequencies: 
low-frequency p modes have their upper turning point rel- 
atively deep in the star and are almost unaffected by the 
treatment of the surface layers, whereas modes of higher 
frequency have turning points closer to the surface and 
the error in their frequencies grows progressively bigger. 
In these circumstances, it is preferable not to calibrate to 
the frequency data themselves but rather to data combi- 
nations chosen to be relatively insensitive to the known 
deficiency in the forward modelling. In the near-surface 
layers, the eigenfunctions of solar p modes of low or inter- 
mediate degree / are essentially independent oil, so the er- 
ror in the frequencies introduced by incorrectly modelling 
this region is just a function of frequency, scaled by the 
inverse of the mode inertia. This suggests that one should 
calibrate frequencies of solar-type stars to data combina- 
tions chosen to be insensitive to such an error. 

Asymptotic analysis suggests other data combinations 
which can be used for model calibration and which are 
more discriminating than the raw frequency data. Mon- 
teiro et al. (2001) consider in some detail the use of the 
so-called large and small separations: 

A„i = i — f„i , Snl = I'nl — l^n-ll+2 (1) 

respectively. As Monteiro et al. demonstrate, the results of 
such a calibration, even for such global characteristics as 




Figure 1. Model calibration selects a model Mo in the 
space of models allowed by the calibration procedure: sub- 
sequent inversion using A4o as reference determines a 
(hopefully improved) model Ai, not necessarily in the orig- 
inal space. If some other calibration selected a different 
model M'q, one hopes that the inversion procedure would 
bring one to the same or similar final model M . 

the mass and age of the star, depend on the other physics 
assumed. 

3. Linearized inversion 

Model calibration is just one approach to the inverse prob- 
lem of inferring the stellar structure, which we may indi- 
cate schematically as X, from the frequency data v. How 
else may the nonlinear dependence i> = v{X) be inverted? 
A way which then allows application of a variety of tech- 
niques is to linearize about a reference model M.^: inver- 
sion takes as input 5v, the difference between the observed 
data and the corresponding values that model Mq pre- 
dicts, and produces as output 5X, the estimated difference 
in structure between the star and A^o- Since the structure 
of M.(j is known, the structure of the star can then be re- 
constructed. Of course, that is a naive hope, because of 
the inherent nonuniqueness discussed earlier, quite apart 
from issues of data errors. But one may at least produce a 
refinement on the initial model and indeed this new model 
can then be used as reference for a subsequent inversion 
as the next step of an iterative approach. Depending on 
the technique adopted, it is helpful and may be essential 
to have a reasonable starting guess in the form of the ref- 
erence model A^o- Model calibration using the large and 
small separations is a reasonable way to find such a model 
for solar- type stars. 
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Model calibration produces a solution A4q that is in 
the span (suitably defined) of the calibration set of models, 
illustrated schematically as a surface in Fig. |^. Inversion 
of the kind just described may then be used to proceed 
from Aio to a new model Ai which may be 'close' to A^o 
but outside the span of the original calibration models. 
One hopes that the final model depends only weakly on 
the initial model, so that if by some other calibration (as- 
suming different physics when calibrating the large and 
small separations, for example) one produces some other 
model A4'q, then the inversion step takes one close once 
more to the same final model A4. 

One can imagine making all structural quantities in 
both the reference model and the target star dimension- 
less by taking out appropriate factors of the gravitational 
constant G and the stars' masses M and radii R. Then 
the frequency differences SuJni between the same mode of 
order n and degree I in the two stars can be related to 
the differences in dimensionless structural quantities by 
an equation such as 



= 6 



<:^^(x)^(x)dx 



F{uJnl) 



(2) 



/o Em 
where here the choice has been made to express the struc- 
tural differences in terms of u, the ratio of pressure p to 
density p, and Y, the helium abundance by mass: for dis- 
cussion of this choice, see Basu et al. (2001). Here, 6u is 
the difference in dimensionless u between the two stars, 
the differences being evaluated at fixed fractional radius 
X in the star. The first term on the right-hand side is 
a constant and just reflects the {GM/ R^Y^^ homologous 
dependence of the frequencies. The final term is some func- 
tion of frequency, divided by mode inertia E„i, which ab- 
sorbs uncertainties from the near-surface layers. For some 
details of the construction of the kernel functions K^J^y i^) 

and iir^"y (x), which are known functions derivable from 
the reference model, see the Appendix. In an inversion, the 
left-hand side will be known; all terms on the right-hand 
side are to be inferred, including the difference in M/R^ 
between the reference model and target star. 

4. Structural and frequency differences 

To motivate the rest of the paper we consider first the 
structural differences between a few stellar models. These 
will already indicate that rather small uncertainties in 
global parameters of stars, e.g. mass or age, can lead to 
large uncertainties in their structure. Of course the posi- 
tive side of that is that the stellar structure is sensitive to 
those parameters and so, if the observable mode frequen- 
cies are in turn sensitive to those aspects of the structure 
then we may have some hope of using the observations to 
constrain e.g. the mass and age of the star rather precisely. 
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Figure 2. Relative differences in pressure p ( solid curve ) 
and density p (dashed curve), at fixed fractional radius, 
between ZAMS .stars of mass I.OMq and I.IMq (in the 
sense I.IMq minus I.OMq). (a) The quantities have not 
been homology-scaled before taking differences, (b) Quan- 
tities have been homology-scaled before differencing. 

Panel (a) of Figure || shows the relative differences in 
pressure p and density p, at fixed fractional radius, be- 
tween two ZAMS stars, of masses I.OMq and I.IMq. We 
note that even for two stars with rather similar masses 
the differences are large, of order unity. Panel (b) shows 
the corresponding differences after the homology scaling 
has been taken out: this scaling is assumed taken out in 
the formulation presented in eq. (||). The effect is essen- 
tially to shift the two curves by a constant: although the 
differences are smaller, they are still large. These changes 
arise from nonhomologous differences in the surface layers 
which change the entropy of the convection zone. Indeed, 
writings = Kp^^^^"- in the stars' convective envelopes (in 
this context only, n denotes polytropic index), and noting 
that sound speed there is essentially determined by surface 
gravity, one finds that in that region 

Ship ~ dlnp 2± — n^lnA'. (3) 

Possibly of more direct relevance for inversion of eq. (||) 
is the difference in adiabatic sound speed c (c^ — jiu, 
where 71 is the first adiabatic exponent). As Fig. ^ shows, 
the relative differences in are smaller than the differ- 
ences in pressure and density, but still quite substantial for 
stars that differ in mass by as little 10 per cent. Homol- 
ogy scaling has little effect on the differences in this case. 
For a reasonable range of masses, these differences scale 
linearly with the mass difference, so for example the cor- 
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Figure 3. Relative differences in sound speed squared (c?) 
at fixed fractional radius, between ZAMS stars of mass 
l.OM© and I.IM© (in the sense I.IM© minus I.QMq): 
the two curves show unsealed (solid) and homology-scaled 
(dashed) differences. 

responding scaled differences between l.OSM© and l.OM© 
ZAMS stars are half those illustrated here, to a very good 
approximation. 

5. Inversion methods 

A number of different linearized inversion methods have 
been developed in helioseismology, geophysics and diverse 
other areas of inversion applications. Two flexible 
approaches are Optimally Localized Averages (OLA) and 
Regularized Least Squares (see e.g. Christensen-Dalsgaard 
et al. 1990 for a description). The OLA method explicitly 
constructs a linear combination of kernels that is local- 
ized at some location in the star and is small elsewhere: 
the corresponding linear combination of relative frequency 
differences is then a measure of the localized average of 
the structural differences. The application of SOLA in the 
asteroseismology of solar-typc stars is discussed and illus- 
trated by Basu et al. (2001). We emphasize here that OLA 
reveals the true extent to which the seismic data alone can 
resolve the aspect of the stellar interior under study. Meth- 
ods which on specific classes of problems appear to have 
superior performance to OLA in resolving aspects of the 
stellar interior are introducing nonseismic information or 
assumptions in addition to the frequency data. 

The form of least-squares method used most exten- 
sively in helioseismology is regularized least-squares: the 
idea is to represent the solution with a set of basis func- 
tions more finely than can be resolved by the data and 
with ideally no bias about the form of the solution built 
into the basis; but then to minimize the sum of chi-squared 
fit to the data and a penalty term which is large if the 
solution has undesirable characteristics. The most used 



penalty function is the integral over the star of the squared 
second derivative of the function under study with respect 
to radius. 

A different way of regularizing the least-squares solu- 
tion, without introducing a penalty term, is to choose a 
drastically smaller set of basis functions. This alternative, 
which we do not claim is intrinsically superior, may give 
apparently better results from few data if the basis func- 
tions are chosen with appropriate intuition or good for- 
tune. The reason is that one can introduce a huge amount 
of prejudice into the solution by forcing it to have a form 
determined by the basis functions. Such a basis could, for 
example, force the buoyancy frequency to be zero in a 
convective core and permit a single discontinuity at the 
core boundary but not elsewhere: such assumptions may 
be reasonable, but it should be realized that they are ad- 
ditional to the seismic data. If the star actually had a 
second discontinuity, or a more gradual variation at the 
core boundary, such an inversion would not generally re- 
veal those features. A helioseismic example is in finding 
the location of the base of the convective envelope, where 
remarkable precision (O.OOli?© taking into account uncer- 
tainties in abundance profiles, O.OOO2if!0 if the only un- 
certainty comes from data noise) has been claimed (Basu 
& Antia 1997): this is credible only insofar as the base of 
the convection zone has precisely the form assumed in the 
inversion, because the true resolution at that location is 
much coarser than that. The true resolution (e.g. verti- 
cally) is essentially limited by the reciprocal of the largest 
vertical wavenumber of the eigenfunctions corresponding 
to the available data (Thompson 1993). 

6. An experiment with specific basis functions 

As a simple illustration of the apparent ability of least- 
squares inversion with a limited basis to infer structure 
even where the mode set has little resolving power, we take 
a basis of five functions chosen for algebraic convenience. 
They are approximately (but only approximately) able to 
represent the difference in u between ZAMS models of 
solar-like stars. The basis represents a function that is 
quadratic in a; = r/i? for a; < 0.1 with zero derivative at 
a; = 0, piecewise constant between radii x = 0.1, 0.3 and 
xi, quadratic between xi and X2 and zero for x > X2- The 
values of xi and X2, which are initially set to 0.65 and 0.75 
respectively, are adjusted by hand to find a minimum of 
the chi-squared fit. 

In the first of two illustrative applications, we con- 
sider the inversion of frequency data (65 modes: I = 0, 
n = 14 - 32; / = 1, n = 13 - 29, / = 2, n = 15 - 30, 
Z = 3, n = 16 - 28) from a 1.O5M0 ZAMS star. In a real 
application we would first calibrate the star by comput- 
ing large and small separations from the data and using 
the approach of Monteiro et al. (2001) to arrive at a ref- 
erence model; in fact we simply took a reference model 
which was a l.OM© ZAMS star. Before discussing the 
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Figure4- Relative frequency differences (v = a;/27r is 
cyclic frequency) between ZAMS stars of masses l.OAfQ 
and 1.05Mq (in the sense I.OSMq minus I.OMq). Lines 
join modes of the same degree: I = (continuous line), 
I — 1 (dotted), I — 2 (dashed) and I ~ 3 (dot-dashed). 
Crosses indicate the actual values of the frequency differ- 
ences, diamonds indicate the values produced by a least- 
squares fit of expression using the basis functions dis- 
cussed in Section 6. . 

results of the inversion, we first consider the data, i.e. 
the relative frequency differences, shown in Fig. ^. (For 
clarity, no noise has been added to the data, though it 
was added for the inversions.) The dominant trends are 
that the values are negative, around —0.066, because of 
the difference in M/R^ between the two stars (in fact, 
S\n{{GM/R^)^) = -0.0668 for these two stellar mod- 
els); and there is a roughly linear trend with frequency, 
in this frequency range, coming from near-surface differ- 
ences. These two contributions can be estimated and re- 
moved by fitting an expression of the form (^: for illus- 
tration, we show in Fig. |^ what would remain: this is the 
signal from the interior, which contains the information 
that the inversion will use to infer conditions inside the 
star. The impression now is that the data contain a signal 
which has some oscillatory component (from the rather 
abrupt change in the structural differences at the base 
of the convective envelope) but is otherwise only weakly 
a function of frequency and increases with decreasing Z, 
indicating that the more deeply penetrating (i.e. lower-Z 
modes) sense 6u/u increasing in the deep interior as one 
gets closer to the centre of the star. 

These features are indeed revealed by the inversion 
(Fig. ^), which compares the exact M-differences with the 
least-squares solutions for noise-free data and for data 
with Gaussian noise with zero mean and uniform stan- 
dard deviations a = 0.1 /iHz, 0.2 fiRz, and 0.3 /iHz. The 
noise realization in the two panels is different, but within 
each panel the noise differs from case to case only by a 
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Figure 5. Residual frequency differences after the fitted 
constant term 6\n{GM/R'^)'^ and a two-term fit for the 
surface-like contribution F{uj)/Eni have been removed. 
Line styles have the same meaning as in Fig. Q. 

multiplicative scaling. It can be seen that for low noise 
levels this very small basis enables the differences to be 
recovered rather well, including the structure beneath the 
convection zone and the variation of Su in the core. For 
larger noise levels the artificiality of the basis functions 
becomes more apparent; also the solution for the higher 
noise levels is rather different for the two noise realiza- 
tions, which gives some indication of the uncertainty even 
in this highly constrained solution. 

A second example is an application to the same mode 
set but with data from model 5*1 of Monteiro et al. (2001) 
and Gaussian data errors with a = 0.1 /iHz. This star 
is slightly more massive and more evolved than the Sun. 
Again, we omitted the calibration step and inverted rela- 
tive to a reference model of the present-age Sun. The re- 
sults are shown in Fig. |7[ Again the qualitative behaviour 
of the differences is recovered reasonably, including the 
downturn in the core. The discrepancies are perhaps partly 
attributable to the fact that the basis functions are not so 
well suited for representing this case as the ZAMS case. 

7. Discussion 

It is a prejudice of some stellar astrophysicists (it was 
indeed expressed a few times in Cordoba during the work- 
shop) that helioseismic experience with the Sun provides 
a poor example when it comes to asteroseismic inversion. 
But the similarities between the two applications are much 
more significant than their differences. An extremely im- 
portant aspect of inversion in any context is to assess what 
the data really tell you and what information is being in- 
troduced by other assumptions or constraints. OLA ker- 
nels (Basu et al. 2001) indicate what resolution can really 
be achieved without additional assumptions. But other 
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Figure 6. Inversion results for the relative difference in 
homology- scaled u between a l.OSM© ZAMS star and a 
I.QMq ZAMS star. The solid, curve shows the exact dif- 
ference du/u, and the short-dashed curve the solution from 
a least-squares inversion of noise-free data. The other 
curves show inversions of noisy data containing indepen- 
dent Gaussian errors with zero mean and uniform stan- 
dard deviation u = 0.1 nHz (dot-dashed), a = 0.2 ^Hz 
(triple dot-dashed) and a = O.i^Hz (long- dashed). The 
upper and lower panel show two different realizations of 
the noise; within each panel the noise from one curve to 
another has just been scaled by a multiplicative factor. 



approaches can advance our knowledge by allowing intro- 
duction of reasonable prejudices: e.g., looking for signa- 
tures of sharp features, or introducing specific basis func- 
tions. The best apparent results are likely to be achieved if 
those functions are physically motivated, because the so- 
lution will accord with onr physical intuition (prcijudicc;). 
This may of course be dangerous. Our example of a highly 
constrained least-squares inversion illustrates that quali- 
tatively reasonable results can be obtained throughout a 
star by introducing assumptions about the form of the so- 



Figure 7. Inversion results for the relative difference in 

homology- scaled u between model Si of Monteiro et al. 
(2001) and a present-day standard solar model (in the 
sense S\ minus solar model). The solid curve shows the 
exact differences, the dashed curve shows the solution of 
the least-squares fit to the data (which contained Gaussian 
noise with uniform standard deviation 0.1 fiHz). 



lution, even in regions where in fact the mode set provides 
little or no localized information (cf. Basu et al. 2001). One 
may have other grounds on which to believe that such a 
solution is plausible, but on the basis of the data alone it 
should be viewed sceptically. 

Model calibration is a useful tool in its own right and 
for obtaining possible starting models for linearized astcr- 
oseismic inversions. Carefully used, model calibration al- 
lows one to build in some prejudice; and the combination 
of calibration and inversion extends the space of solutions 
that one explores. 

Calibrating on large and small separations can be ef- 
fective (Monteiro ct al. 2001), always assuming that one 
has not neglected some important aspect of the physics: 
in that regard, the effects of rotation and magnetic fields 
need to be borne in mind (sec Dziembowski & Goupil 
1998), in solar-type stars as in many other pulsating stars. 

Finally we note that mode identification may be a 
problem, even in solar-type stars. The asymptotic pattern 
of high-order p-mode frequencies may allow / to be deter- 
mined, but there may be some uncertainty in n. Such an 
uncertainty can be allowed for in the inversion by adding 
to the right-hand side an extra term 

G{uJnl)/u}nl (4) 

where G is some function of frequency. We have verified 
that, with our ZAMS example, even restricting G to be 
a constant function removes quite satisfactorily the effect 
of a misidentification of n. More generally, G could ju- 
diciously be chosen to reflect the variation of the large 
separation with frequency. 
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Appendix 

The derivation of kernels relating the linearized differ- 
ences in structure to the differences in frequency (cf. eq. |^) 
has been discussed by, e.g., Gough & Thompson (1991), 
Gough (1993), and Kosovichev (1999). As all those au- 
thors show, the kernels for either of the pairs of variables 
(c^,p) or (71, p) are quite straightforward to derive from 
the equations of linear adiabatic oscillations together with 
the linearized equation of hydrostatic support. Obtaining 
kernels for various other pairs, including the pair (u, Y) 
used in this paper, can be accomplished by first obtaining 
kernels for one of the other two pairs and then using the 
following piece of manipulation. It is sufficiently ubiqui- 
tous (occurring often when one wishes to transform from 
kernels for a pair including p to some other variable pair) 
that we write it rather generally. 



Let ip{r) be a solution of 
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(5) 



where prime denotes differentiation with respect to r (or x 
if everything - including r - is expressed in dimensionless 
variables), for a given function F{r) (so -0 is a functional 
of F), with boundary conditions '(/'(O) = and ip{R) = 0. 
Then, provided 5m{Q) — and 5m{R) = 0, where mir) 
is the mass interior to radius r, and letting (...) denote 
integration from r = to r = i? (a; = to x = 1), 



P / \ u 



(6) 



Note that 6m{0) = holds in general; and in our present 
application, we scale all structural quantities by G, M and 
R to make them dimensionless, so 6m{R) = is forced to 
be true. These two conditions also mean that (^Airr'^dp) is 
zero, so any multiple of 47rr^p may be added to a kernel 
multiplying 5p/ p. Note that such an additional contribu- 
tion to F{r) makes no change to the right-hand side of 
eq. (|^) and hence the contribution to ij] from such an ad- 
dition is zero. 

Obtaining kernels for Y additionally requires an as- 
sumption about the equation of state through 71, since 
the oscillations do not know directly about the chemical 
abundances, hi the following we write 

/ain7i\ /^<91n7i\ 



Then for convenience we record the following transforma- 
tions: 



(c^p) (m,7i) 



(8) 



with F = Kp^2] 
(71, p) ^ iu,Y) 



Ku,Y =l,pK^^^p-p(^] (9) 



with F = (7,p +7,p )K^^^p + Kp^^^; and 
(m,7i) -> {u,Y) 

Ky.u = 1,Y 



Ku,Y = l,pK^^^u + Ku^y I p ) ^-^^^ 



with F = (7,p -1-7, p )i^^i,„ 



